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Abstract 

We discuss the inclusion of running coupling effects in perturbative small x evolution 
equations. We show that a running coupling BFKL-like x-evolution equation is fully 
compatible, up to higher twist corrections, with the standard factorized perturbative evo- 
lution of parton distributions. We then use this result, combined with the well-known Airy 
asymptotics, to prove that the oscillations which are present in the running-coupling BFKL 
solution do not affect the associated splitting functions, which instead remain smooth in 
the small x limit. This allows us to give a prescription to include running-coupling cor- 
rections in the small-x resummation of scaling violations. We show that these corrections 
are small in the HERA region. 
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1. Introduction 

The theory of scaling violations of deep inelastic structure functions at small x has 
recently attracted considerable interest, prompted by the extensive results obtained by 
experiments at HERA [1,2]. New effects beyond the low-order perturbative approxima- 
tion [3-5] to anomalous dimensions or splitting functions should become important at 
small x. The BFKL approach [6-10] provides in principle a tool for the determination of 
the small x improvements of the anomalous dimensions [11,12]. However, the data are in 
good agreement with a standard next-to-leading order perturbative treatment of scaling 
violations [13,14,1,2], while naive attempts [15] to incorporate small x logarithms fail com- 
pletely [16,17]. Important progress has been made recently in understanding the empirical 
softness of the BFKL resummation [18-21]. Specifically, it is now understood that in order 
to avoid instabilities [22,23], the inclusion of small x contributions must be combined with 
a resummation of the collinear singularities of standard perturbative evolution. 

An important aspect of the problem has to do with the implementation of running cou- 
pling effects. Indeed, it is a priori not obvious that the 'duality' relations [11,24,19] which 

connect BFKL evolution in £ = In - and renormalization-group evolution in t = In ^ 
remain valid when the running of the coupling is taken into account. These duality re- 
lations play a crucial role [19,20] in the desired simultaneous resummation of collinear 
singularities and small x logs, and imply that in small-x evolution factorization is pre- 
served, namely, anomalous dimensions remain independent of the boundary conditions on 
the parton density. 

In our previous work [23,19,20] we have included running coupling effects pertur- 
batively up to NLLx. We have shown that duality and factorization continue to hold, 
provided the anomalous dimension is supplemented by a term induced by the running of 
the coupling [25,23]. This term however corresponds to a singular contribution to the 
BFKL kernel, or equivalently, a perturbatively unstable contribution to the anomalous 
dimension. In MS factorization this singular contribution is shifted to the quark sector [7] . 
Even though its effects are in practice rather small in the HERA region [20] , this singular- 
ity signals a failure of the NLLx treatment of the running of the coupling in the asymptotic 
small x limit. 

In this paper, we address in detail the problems related to the inclusion of running 
coupling effects in small x evolution equations. After reviewing the way duality is affected 
when the running of the coupling is included perturbatively, we consider the modified 
running-coupling BFKL x-evolution equation and its solution. We discuss the properties 
of the corresponding anomalous dimension and prove factorization and duality to all per- 
turbative orders, up to higher twist corrections. 1 We then study in detail the running 
coupling contribution to anomalous dimensions in the asymptotic small-x limit, which may 
be determined exactly in terms of Airy functions, by taking a quadratic approximation of 
the BFKL kernel near its minimum [27]. In particular, we prove that the unphysical oscil- 
latory behaviour which has been shown [21,27-32] to characterize solutions to the running- 
coupling BFKL equations at small x does not affect the small x asymptotics of splitting 

1 Arguments for factorization have already been presented elsewhere [21,26], but they are 
partly based on model-dependent assumptions. 
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functions. Therefore, the perturbative approach to scaling violations need not break down 
in the small x limit. We then show that this quadratic approximation is sufficient to de- 
termine the leading x — > asymptotics, since only the behaviour of the kernel near the 
minimum is relevant in this limit. Finally, we discuss how this asymptotic behaviour can 
be systematically matched to the small x expansion of anomalous dimensions, and thus 
used to resum small x running coupling terms. We show that the resummed result is free 
of unphysical singularities, perturbatively stable, and smooth in the small x limit, and 
in the HERA region it only leads to a small correction to our previous phenomenological 
results. 



2. Perturbative Duality 

We start by reviewing the way duality of perturbative evolution is proven, first at 
fixed coupling and then perturbatively in the running coupling case. We then show how 
the perturbative expansion in the running coupling case is beset by unresummed singular 
contributions. The behaviour of structure functions at small x is dominated by the large 
eigenvalue of evolution of the singlet parton component. Thus we consider the singlet 
parton density 

t) = x[g{x, Q 2 ) + k q ® q{x, Q 2 )], (2.1) 

where £ = logl/x, t = log (J 2 / ^ 2 , g(x,Q 2 ) and q(x, Q 2 ) are the gluon and singlet quark 
parton densities, respectively, and k q is such that, for each moment 

/■OO 

G(N,t)= / dfe-^Gfcf), (2.2) 
Jo 

the associated anomalous dimension 7 (ct s (£), A) corresponds to the largest eigenvalue in 
the singlet sector. The generalization to the full two-by-two matrix of anomalous dimen- 
sions is discussed in detail in ref. [20] . 

At large t and fixed £ the evolution equation in A-moment space is then 

j t G(N,t)=^(a s (t),N)G(N,t), (2.3) 

where a s (t) is the running coupling. The anomalous dimension is completely known at 
one- and two-loop level: 

7 (a s , A) = a s7o (A) + a 2 7i (A) + . . . . (2.4) 
The corresponding splitting function is related by a Mellin transform to 7 (ct! s , A): 

1 (a s ,N)= [ dxx N P{a s ,x) (2.5) 
Jo 

Small x for the splitting function corresponds to small N for the anomalous dimension: 
more precisely P ~ l/x(log(l/x)) n corresponds to 7 ~ n\/N n+1 . Even assuming that a 
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leading twist description of scaling violations is still valid in some range of small x, as 
soon as x is small enough that a s £ ~ 1, with £ = \ogl/x, all terms of order (a s /x)(a s £) n 
(LLx) and a s (a s /x)(a s ^) n (NLLx) which are present in the splitting functions must be 
considered in order to achieve an accuracy up to terms of order a 2 .(a s /x)(a s £) n (NNLLx). 
As is well known, these terms can be derived from the knowledge of the kernel 

X (a s , M) = a 8Xo (M) + a 2 sX i(M) + ■■■ ■ (2.6) 

of the BFKL ^-evolution equation 

^G(£,M) = x(a s ,M)G(£,M), (2.7) 
which is satisfied at large £ by the inverse Mellin transform of the parton distribution 

/oo 
dte- Mt G(^t). (2.8) 
-oo 

This derivation was originally performed [11] at LLx by assuming the common validity of 
eq. (2.7) and eq. (2.3) in the region where Q 2 and £ are both large. However, it was more 
recently realized [24,19,20] that the solution of eq. (2.7) coincides generally with that of 
eq. (2.3), up to higher twist corrections, provided only that the kernel of the former is 
related to that of the latter by a 'duality' relation, and boundary conditions are suitably 
matched. This implies that the domains of validity of these two equations are in fact the 
same in perturbation theory, up to power-suppressed corrections. 

The derivation of duality is simplest when the coupling does not run, in which case 
the relation between the kernels of the two equations is 

X [a 8 , 1 (a 8 ,N)] = N. (2.9) 

This result is immediately found by writing both eq. (2.3) and (2.7) and their solutions in 
terms of the inverse Mellin transform of the parton distribution 

/•OO />oo 

G(N,M)= d^e~ N ^G(^M)= dt e~ Mt G(N,t). (2.10) 

JO J -oo 

The inverse Mellin of the solution to eq. (2.7) is 

where Fq(M) is an iV-independent boundary condition. The large t behaviour of G(N,t) 
is determined by the rightmost singularity of G(N,M) (2.11) in the M-plane, while the 
contributions of additional singularities further to the left are suppressed by powers of Q 2 . 
In particular, perturbative singularities are given by solutions of the equation x( a s, M) = 
N, while singularities of the boundary condition F (M) are nonperturbative. 
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Expanding the denominator of the solution (2.11) about its rightmost perturbative 
singularity, the solution eq. (2.11) is thus seen to coincide, up to higher twist and nonper- 
turbative terms, with the inverse Mellin of the standard solution to the renormalization 
group eq. (2.3), namely 

G(N,M) = (2.12) 

1 ' M-j(N,a a ) 1 ' 

provided the duality equation (2.9) holds, and the t-independent boundary condition 
F (N) to eq. (2.3) is related to F (M) by 

F0(N) ~ -X'h(*s,N))- (2 ' 13) 

This establishes the perturbative equivalence (duality) of eqs. (2.7) and (2.3), up to higher 
twist corrections. Using the perturbative expansion of x (2-6) in the duality relation (2.9), 
we then find that %o determines n f s (a s /N), while a s xi leads to a s "f ss (a s /N), where 
the combination r y s (a s /N) + a s "f ss (a s /N) correspond respectively to all terms of order 
(a s /x)(a s {;) n and a s (a s /x)(a s ^) n in the splitting functions. 

When one goes beyond LLx, i.e. beyond the leading-order approximation for x, the 
running of the coupling cannot be neglected, and this derivation must be re-examined. 
Indeed, in M space the usual running coupling a s (t) becomes a differential operator: 
taking only the one-loop beta function into account 

Q, 

= - ' d +••-, (2-14) 

where (3q is the first coefficient of the /3-function (so (3 = —f3oa 2 s + " 0; w ^ n the obvious 
generalization to higher orders. Hence, assuming the coupling to run in the usual way with 
<5 2 , the ^-evolution equation eq. (2.7) becomes [27] 

^G(£, M) = X (a s , M)G(£, M), (2.15) 

where the derivatives with respect to M act on everything to the right, and x ma y be 
expanded as in eq. (2.6) keeping the powers of a s on the left. Clearly, an operator-valued 
evolution kernel is a potential threat to factorization. Furthermore, one may ask whether 
eq. (2.15) is indeed the proper form of the running-coupling BFKL equation. Here we will 
show that, besides being intuitively appealing, the evolution eq. (2.15) has the property 
that its solutions are still the same as those of a renormalization-group equation eq. (2.3), 
with a suitably matched kernel, to all perturbative orders and up to power corrections. This 
implies consistency of eq. (2.15) with standard perturbative evolution, and in particular 
that factorization is preserved. Therefore, eq. (2.15) can be viewed as an alternative 
representation of the standard renormalization-group equation. 

It is clear from eq. (2.15) that running coupling effects begin at NLLx. In fact, to 
N fe LLx it is sufficient to retain the first k terms in the expansion of a s (2.14) in powers of 
a s ; the ensuing equation can then be solved perturbatively. To NLLx one finds [23] that 
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this solution is again the same as that of a dual t-evolution equation (2.3), provided the 
fixed-coupling duality relation eq. (2.9) is modified by letting a s — > a 8 (t), and then by 
adding to 7 SS an extra term [25] A7 SS proportional to fto: 

A^ 8S ( w )--Po 2x , 2{is) • (2.16) 

Equivalently, the duality relation eq. (2.9) can be formally preserved, provided that a s — > 
a s (t) and the function x used in it is no longer identified with the BFKL kernel, but rather 
given by an 'effective' x function 

Xeff(« s ,M) = x(a s ,M) + A X (a s ,M), 

(2.17) 

A x (a s ,M) = a 2 s A X i(M) + ..., 

where x(a s , M) is obtained letting a s — » a s in the kernel of eq. (2.15), 2 and the correction 
term A% to NLLx is given by 

Axi(M ) = A^^ge. (2.18) 

Also, the matching of boundary conditions is now given by a more complicated relation 
than eq. (2.13), involving derivatives of Fq and xo'i however, it remains independent of 
t, and only involves N and a s (at the initial scale). We have further verified by explicit 
computation that if the perturbative solution of eq. (2.15) with the running coupling 
eq. (2.14) is pursued up to N 2 LLx these results remain valid: the boundary condition 
remains t-independent, and duality can be preserved, provided only a further Ax2 term 
is included in Xe« eq. (2.17). 

The problem with the perturbative approach is that the correction terms which must 
be included in the effective x e Q- (2-17) have an unphysical singularity: xo(M) has a 
minimum at M = -|, so the denominator of Axi(M) vanishes, resulting in a simple pole 
in the NLLx correction Axi eq. (2.18) at M = |. The NNLLx correction A% 2 turns out 
to have a fourth-order pole, and in fact at each extra order three extra powers of (xq) _1 
appear. As explained in ref. [23], this leads to a perturbative instability: as a consequence 
of the singularity, the splitting function AP SS eq. (2.5) associated with the anomalous 
dimension A7 SS behaves as 

AP ss (a s ,Q 2 
p / n ~ • ( 2 - 19 ) 

In practice, the coefficient of this rise turns out to be small enough that its effects are 
negligible in the HERA region [20]. However, the situation remains unsatisfactory from a 
theoretical point of view: it would be better to extend the perturbative proof of factor- 
ization and duality to the case when the running of the coupling is included to all orders, 
and then sum up these perturbatively unstable terms. In the sequel, we will do this, and 
show that the M = \ singularity and the corresponding ones that appear at higher orders 
in a s (3o are an artifact of the expansion, are not present in the all-order solution, and can 
thus be resummed. 



2 Notice that, beyond leading order, this is not the same as the fixed-coupling kernel of 
eq. (2.7), which refers to the unphysical case of a theory where 0o = identically. 



5 



3. Running Coupling Duality and Factorization 



In this section, we prove to all orders the perturbative duality and factorization in- 
troduced in the previous section, and we discuss some general properties of the all-order 
solution to the running coupling BFKL equation, specifically its apparently unphysical 
behaviour in the £ — > limit. We start from the running coupling £ evolution equa- 
tion (2.15). We consider first the case in which the running is included at one loop accord- 
ing to eq. (2.14), and the kernel x is linear in the operator a s : 



x{a s ,M) = a s ip{a s ,M), 



(3.1) 



where ip is a function of M and the fixed coupling a s . This is of course the case when the 
leading order form of x is considered, in which case <p = xo(M), but also more generally 
whenever terms of higher order in the fixed coupling a s are included in ip, as in the case of 
the resummed leading order kernel xo previously discussed by us [19,20]. In the remainder 
of this section we will omit the explicit dependence of ^ on « s , since it will not play any 
role in the ensuing proof of factorization and duality. The generalization of the proof to 
higher orders in the running coupling a s (t) will be discussed at the end of the section. In 
sect. 5 we will then show that, in order to perform a running coupling resummation to 
NLLx, it is sufficient to consider a kernel x of the form of eq. (3.1), and we will also discuss 
higher order generalizations. 

With the kernel (3.1), after taking a second Mellin transform (2.2), eq. (2.15) becomes 
the differential equation 



(l-P a s ^)NG(N,M) + F(M) 



a s cp{M)G{N, M) 



(3.2) 



The function F(M) is the boundary condition, obtained by acting with the operator (1 — 



/3 a t 



dM- 



on the boundary condition — Fq(M) of eq. (2.12). The perturbative solution [23] 



to eq. (3.2) discussed in the previous section is an expansion in a s at fixed a s /N and it can 
be built up by an iterative procedure: at order n the order n — 1 expression for G(N, M) 
is inserted in the right-hand side. Clearly, this iterative solution is a linear functional of 
F(M). 

However, it is also easy to determine the general solution to the differential equa- 
tion [27,28]: 



G(N,M) = H(N,M) + [ M dM> H ^ M \ 
K ' ) y } J mo H(N,M>) [3 a s N' 



(3.3) 



where H(N, M) is the solution of the homogeneous equation 



H(N, M) = H(N, M ) exp 



M -M ( 



o 



foot i 



(3 N 



M 



dM'ip{M') 



M 



(3.4) 



with an arbitrary initial condition H{N, M ). The inhomogeneous term can be written as 



pM 

G(N, M) = / ( 

J M 



dM' exp 



M-M' 

(3 a s 



M 



dM"ip(M") 



M' 



F(M') 



(3.5) 
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and thus depends linearly on F(M) while being independent of H(N,M ). 

Let us now see how the iterative solution can be recovered from the general solution, by 
expanding it in a s at fixed a s /N: the iterative solution corresponds to the inhomogeneous 
term in eq. (3.5), while the homogeneous solution in eq. (3.4) vanishes faster than any 
fixed perturbative order in (3q. To see how this works, change the integration variable in 
eq. (3.3) from M' to 

i r M 

y(M, M') = — / [N - a sV {M")]dM" (3.6) 

(J a s N J M , 

The inhomogeneous solution (3.5) then becomes 

G(N, M)= f 
Jo 



v(m,m 0) F (M'(y,M)) 

[N-aMM>(y,M))f y > 



where M' is implicitly defined as a function M' = M'(y,M) by eq. (3.6). The expansion 
of M'(y, M) in powers of a s is 

M '= M+ w^m v+oia2 - ) - (3 - 8) 

Using this expansion to do the integral, and noting that y — > oo in the limit a s — > 0, one 
finds 

G(NM )= F{M) I asPo d F{M) I (39) 
{ ' ' [N — a s cp(M)] [N — a s cp(M)] dM [N — a s cp(M)] ' ^> 

which is the iterative solution of ref. [23], discussed in the previous section. Furthermore, 
the homogeneous solution (3.4) is manifestly proportional to e~ y and therefore vanishes to 
all orders in the perturbative expansion. Thus in the following we drop the homogeneous 
term and study only the inhomogeneous term (3.5), which reproduces the perturbative 
solution. 

Starting from G(N, M) one can do either of the inverse Mellin transforms, and ob- 
tain G(£,M) or G(N,t). We consider G(N,t) first, in order to prove that it satisfies a 
renormalization-group equation of the usual form, and study the modification induced 
by running-coupling effects on the anomalous dimension and the duality equation which 
relates it to the kernel \- By taking the inverse M-Mellin transform we can write G(N, t) 
in the form 

G(N , t) = exp [Mt + Mj^o _ 1 r Vv<An] (3.10) 

Jc-ioo 2tu L (3 a s (J N J Mq J N 

where 

I(N,M)= r\lM'e,p\^-^ + 1 ±- f M ' dMMM")}^ 1 - (8.11) 
J m 1 Po« s Po^v J Mo J fj a s 
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It is well-known [28,25] that evaluating the M integral by saddle point one recovers a 
running-coupling version of the duality relation eq. (2.9). However, both duality and 
perturbative factorization seem to be spoiled in the process. We now review this derivation, 
and then prove that, contrary to appearances, factorization still holds in the perturbative 
limit, and duality is recovered with a series of corrections eq. (2.17). 

The saddle-point condition for the integrand of eq. (3.10) (assuming I(N, M) to be a 
smooth function of M) is 

fj a s (J N 

Noting that at leading order t + l/(/3oa s ) = [Poa s (t)] , we see that this is equivalent to 

a s (t)<p(M s ) =N. (3.13) 

Identifying M s with the anomalous dimension, 

M s = ls (a s (t)/N), (3.14) 

we see that eq. (3.13) coincides with the duality relation, eq. (2.9), but with the fixed 
coupling replaced by the running coupling [28,25]. Substituting back the saddle condition 
in the exponent of eq. (3.10) we get (differentiating the duality relation with respect to t) 

M s Mo 1 



/ dMV(M') = / dt' ls {a s {t')/N), (3.15) 

a s (t)/Jo a s p Poa s J\ J Mq J to 

where r y s (a s (t')/N) is defined for all a s (t') as the solution of eq. (3.14), and we choose for 
simplicity (and with no loss of generality) a s such that a(to) = a s . Performing the saddle 
integral we thus find 



G(N, t) = y _^ [Mj(t)] exp ^ dt ls (a(t ), N) J + • • • , (3.16) 

where the dots denote corrections to the leading saddle-point approximation. 

Eq. (3.16) seems to violate factorization because of the t-dependence induced by the 
substitution of the saddle condition M = M s (t) in the integral I(N,M), which depends 
on the boundary condition. Also, it seems to violate duality because the t-dependent 
prefactor 1/y 7 ^ 7 spoils the identification of the anomalous dimensions with 7 S . This is 
puzzling because these violations seem to start at NLLx, whereas the explicit perturbative 
calculation discussed above shows that factorization and duality are respected at NLLx. 
We now address both issues in turn. 

First, let us consider the integral I(N, M s (t)) which appears in eq. (3.16) (with I(N, M) 
defined in eq. (3.11) and M s eq. (3.14)) in the perturbative LLx and leading-log Q 2 limit, 
i.e. as a s — > 0, with a s /N and a s t fixed. We change the integration variable from M' to 
u(Mq, M') defined in eq. (3.6), so we can write this integral in the form 

ryh(a(t)/N)] F(M'(v,Mo)) 



8 



We take (1 — jf<fi) < and M' < M , so that y > 0. As a result, in the perturbative limit, 
nil) ~^ 00 an d the perturbative t dependence disappears to all orders in the expansion. So 
perturbatively I(N, M s (t)) is independent of t to all orders, even though it does depend 
on to because M' = M'(y, Mq). This proves factorization to all perturbative orders: the 
t dependence of G(N,t) is entirely determined by the evolution kernel %, and does not 
depend on the boundary condition 

The t-dependence, however, is not only due to the anomalous dimensions , y s (a s (t)/N) 
which satisfies the running-coupling duality relation eqs. (3. 13), (3. 14), but also to the 
square-root prefactor in eq. (2.16). However, because this additional t-dependence goes 
entirely through the running coupling a s (t) , it has the structure of a NLLx coefficient func- 
tion, and it can thus be reabsorbed by small-x scheme change [23,20,33] in a redefinition 
of the anomalous dimension and the boundary condition: 



V(7,(^#)) 



exp — 



IlnV(7 s (^#)) 



exp 



exp 



t ///q s (t')\ /q s (t')\ 
Po I dt' — — — ~T7rr^ — Oi s {t') 
V 2 (^) 



-^(^r 2 )) 



(3.18) 



dt' A lss (a s (t')/N) 



to 



V(7s(^)) 



where the second equality follows from differentiating the duality relation, and the last 
factor, which depends on the initial scale to only, can be reabsorbed in the boundary 
condition. Identifying <p = xo we see immediately that the resulting addition to 7, A7 SS , 
is the same as that determined perturbatively at NLLx eq. (2.16). 

Furthermore, we can use the saddle-point method to generate an asymptotic expan- 
sion of G(N,t). Noting that the first correction to the gaussian approximation is a cubic 
term in M — M s proportional to ip', and that odd terms in the expansion vanish upon 
integration, this is seen to generate an expansion in powers of 



N(3 



,3 



= a. 



a 



^) h' s (a s (t)/N)f[3 , 



(3.19) 



where 7^ is differentiated with respect to its argument a s (t)/N. It follows that higher 
orders in the expansion give a series of terms which can be absorbed in a series of contri- 
butions to the anomalous dimension of the form 



7rc (a s (t),A0 = 7s (- 



.(*)■ 



+ /3 a s (t)A 7ss (^#) + [/3 a s (t)] 2 A 7sss (^) + 



N 



(*)> 



(3.20) 



where A7 SS , given by eq. (2.16), is due to the leading fluctuations about the saddle, A7 SSS 
comes from the next-to-leading fluctuations and so on. It is also easy to check that the 
boundary condition does not contribute to this additional t dependence, because all its 
higher order derivatives with respect to M vanish in the perturbative limit, for the same 
reason that the homogeneous term vanishes to all perturbative orders. 
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This shows that we can formally preserve the duality relation eq. (2.9), with a s — > 
ct s (t) and an effective Xefi e Q- (2-17), provided at each order in a s we add a new correction 
term to the kernel determined from 7 rc eq. (3.20); the first correction is given by 
eq. (2.18). The order a£ correction is proportional to [<^'] _3fc eq. (3.19), and thus if cp has 
a minimum, the correction has a 3/c-th order pole as M — > M min . Therefore [23], when 
£ — > oo the corresponding splitting function rises by a factor £ 2fc faster than the leading 
order, and the perturbative expansion is unstable. These singular contributions will be 
resummed in sect. 4-5, leading to an anomalous dimension 7 and associated kernel XeS 
which are free of singularities, and thus have stable perturbative expansions. 

Before we turn to this task, let us exploit further the factorization property which 
follows from the fact that the integral I(N, M s (t)) in eq. (3.17) is asymptotically a function 
of N, independent of t to all orders in the perturbative expansion in a s (t). This implies 
that, in the same limit, G(N,t) in eq. (3.16) can be written as 



/•c+ioo 

G(N, t) = / 

J c—ioo 



dM 
2ttT 



Mt 



exp 



M — Mq 

(3 a s 



M 



dM"ip(M") 



Mn 



I(N,t ) 
N 



(3.21) 



Note that this is the same form that one could obtain from G(N, M) given by the homo- 
geneous solution in eq. (3.4) with the identification H(N,Mq) = I(N, to)/N, that is with 
an initial condition which is a linear functional of F{M). 

This proof of factorization and duality can be extended to the case in which the run- 
ning of the coupling is included to higher orders. To this purpose, introduce a perturbative 
expansion of the inverse M-Mellin G(£,t) of the solution to eq. (2.15) in powers of a s (t): 



<?(£, t) = G (£, t) + (£, *) + .... (3.22) 

Then, the double Mellin Gq{N,M) satisfies eq. (3.2) with ip = xo, while G\ satisfies 



NG 1 (N,M) = a s [ X0 (M)G 1 (N,M) + xi(M)G (N,M)} + O(a 2 s ), (3.23) 

where we have made use of the fact that the commutator of a s and %o is of order a 2 s . 
Eq. (3.23) is immediately recognized to have the same form as eq. (3.2), with a boundary 
condition now given by the XiGq term. Hence, its solution has the form eq. (3.3) and 
factorization and duality are preserved. Specifically, G\ is obtained multiplying the factor- 
ized Go eq. (3.21) by a suitable coefficient determined by x\i an d thus it has the structure 
of the standard perturbative solution, but now with the factor of a s which multiplies xi 
promoted to running a s (t). The procedure can then be pursued to higher orders: in fact, 
thanks to the expansion of G in powers of a s (t) , the higher-order differential equation for 
G which is found by simply letting a s — > a s in the higher orders of the expansion of x 
eq. (2.6) is reduced to a linear differential equation for each individual Gi. 

Finally, it is interesting to observe that the inverse iV-Mellin transform of G(N, M) 
can actually be evaluated exactly if the boundary condition I(N, M) is assumed to be 
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a sufficiently smooth function of N, so that it may be effectively considered to be N- 
independent, i.e. assuming that the boundary condition is ^-independent, and all further x 
dependence is generated by perturbative evolution. Indeed, using for definiteness eq. (3.21) 
and setting I(N, to) = Goo one gets [34] 




G(g,M) =Gooexp M " ^ u / I 2W£- / dM'<p(M') I , (3.24) 



M 
M 



where I is the modified Bessel function. Even though it is not possible to determine 
the inverse M-Mellin transform of eq. (3.24) in closed form, the asymptotic behaviour of 
G(£,t) when t is large may be determined by performing the inverse M-Mellin by saddle 
point. The saddle condition is 

1 " (M » (3.25) 




In the limit of large £, the left-hand side of eq. (3.25) vanishes, while on the right-hand 
side Iq/Io ~ 1 so the whole expression vanishes only if ip goes through zero for some 

M (recall that we assume M < M in order for the perturbative limit to exist). Clearly 
if (f is always positive, as is the case for ip which has a minimum with a positive value, 
the saddle equation cannot be satisfied at a real value of M: in such case, the asymptotic 
behaviour of G(£, t) is dominated by a pair of complex saddle points, and therefore displays 
unphysical oscillations. Alternatively, if ip has a vanishing or negative minimum, there is 
always a real saddle. The asymptotic limit is then stable, provided we choose M equal 
to the value of M where x vanishes (otherwise, for negative minimum the square-root in 
eq. (3.24) becomes complex). 

The onset of unphysical behaviour at large enough £ has been interpreted as a sig- 
nal of breakdown of perturbation theory [21,28,30-32], and thus it has been viewed as 
an indication that the perturbative approach must be abandoned in the small x limit. 
However, in the following section we will show that if we are interested in the BFKL ap- 
proach as a tool to improve the theory of scaling violations at small x, then this negative 
conclusion is unjustified. Indeed, even when G(£, t) determined by inverse Mellin from 
eq. (3.24) displays an asymptotic oscillatory behaviour, the associated splitting function 
remains smooth and monotonic, and therefore no oscillations are seen in the factorized 
solutions to leading-twist evolution equations. 



4. Quadratic Kernel and Asymptotic Behaviour 

In order to resum the running-coupling corrections to duality eq. (3.20) we must 
determine the parton distribution G(N,t) by performing the inverse M— Mellin transform 
in eq. (3.10) exactly, rather than in the saddle-point approximation as in the previous 
section. Even though it is nontrivial to do this in general, the integral can be easily 
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determined [27,25] if the kernel cp in eq. (3.10) is quadratic. The ensuing form of the parton 
distribution displays the unphysical oscillations discussed at the end of the previous section. 
In this section we review this determination of G(N, t), and then we use it to determine the 
running-coupling dual anomalous dimension eq. (3.20), and show that the corresponding 
splitting function, and thus the solution to the associated evolution equation, are free of 
oscillations and have a stable asymptotic £ — > oo limit. In the next section we will then 
show that knowledge of the resummed anomalous dimension for a quadratic x is sufficient 
to resum running-coupling singularities in general. 
We start assuming that <p is given by 



<p(M) =c+ -k(M- \] 



(4.1) 



where k and c are constants. In the sequel we will discuss the three cases of positive, 
negative, or vanishing c. While we will assume that n is positive but otherwise generic, for 
explicit numerical examples we will take 



K = — 



2£a „/,"/■ I 



in 



(4.2) 



(where ip(x) is the digamma function and Ca = 3) which corresponds to the curvature 
of the leading-order BFKL kernel [6] xo (2-6) around its minimum. The expression of 
G(N, t) is then given in terms of an Airy function [27]: inserting eq. (4.1) in eq. (3.21) and 
performing the M integral with the choice M = |, we get 



G(N,t) =K(iV)e 1 /' 2/3oa -W]Ai[z(a s (t),iV)] /(i ^ to) . 

where Ai(z) is the Airy function, which satisfies 

Ai"(z) - zM(z) = 0, 

with Ai(0) =3" 2 / 3 /r(2/3), 



(4.3) 



(4.4) 



z{a s {t),N) 



2f3 N\ 



1/3 



K 



1 



a s (t) N 



(4.5) 



and K(N) is a t-independent normalization, 



(4.6) 



Along the positive real axis, the Airy function is a positive definite, monotonically 
decreasing function of its argument, and it behaves asymptotically as 



Ai(z) 



iTT- 1 /^- 1 / 4 



exp(-f^/ 2 )(l + 0(z- 3 / 2 )). 



(4.7) 
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0.075 0. 1 0. 1250 . 150 . 175 0. 2 0. 225 

Figure 1: Contour plot in the complex iV plane for the Mellin inversion integral of 
G(N,t). The crosses indicate the approximate location of the complex saddles, and the 
dotted line is the approximate steepest descent path. 

For large enough iV ', z ~ iV 1 / 3 , so Ai[z(a s , N)] is also a monotonically decreasing function 
of AT, asymptotically damped as Ai[z(a s , AT)] ~ exp—y/N. When c = 0, this monotonic 
behaviour persists for all real positive iV even when iV is small. However, if c > 0, z changes 
sign when iV < a s (t)c, and Ai(z) oscillates along the negative real z axis, so Ai[z(a s , N)] 
oscillates as A" — > 0, because z — > — oo in this limit. If instead c < 0, z — > oo not only as 
A" — > oo but also as A^ — >■ 0, and Ai[z(a a , N)] is smoothly damped in both limits. 

It is now easy to check that the behaviour of G(£, t) obtained by inverse Mellin trans- 
formation from G(N, t) (4.3) has the asymptotic properties discussed in general at the end 
of sect. 3. Specifically, if c (i.e. the value of x (4-1) at its minimum) is positive, then 
G(£,t) displays an unphysical oscillatory behaviour as £ — > oo. This can be seen explicitly 
by evaluating by saddle point the Af-Mellin inversion integral which gives G(N,t). The 
saddle point condition is dominated by the Airy function, and has thus the form 

dM[z(a s ,N)]/dN 
? Ai[z(a.,N)] ' 1 ' ] 

If £ is not too large, the saddle condition is satisfied in the large A^ region where Ai[z(a s , N)] 
is a decreasing function of N. As £ grows, the saddle is drawn towards A^ — >• where 
Ai[z(a s , N)] decreases faster with A^; however, for sufficiently small A^ the derivative 
changes sign and the real saddle is lost. This situation is displayed in fig. 1, where we 
show a contour plot in the complex A^-plane of the real part of the integrand to the A^- 
Mellin inversion, namely Re [£N + ha.G(N, t)] with G(N, t) given by eq. (4.3), I(N, to) = 1, 
and t = 6 (corresponding to a s (t) ~ 0.2), £ = 10, c = 1. Along the real axis the function 
increases, and there is no saddle. However, even when £ — >• oo the saddle condition eq. (4.8) 
does still have a pair of complex solutions, with |AT| — > as £ grows. These saddle points 
can be clearly seen in fig. 1. 

For large £, the position of the complex saddle points can be determined explicitly 
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using the asymptotic expansion eq. (4.7), which is valid when | arg(z)| < 7r: one finds 

N± = e ±l7r/4 (k/0 1/2 , (4.9) 

where k = f^-(^) 1 ^ 2 - The asymptotic large £ behaviour is dominated by the contributions 
from this pair of complex conjugate saddle points, which add to give cosine oscillations: 
explicitly 

t) exp( cos (v^e) . (4.10) 

As £ becomes large, G(£,t) becomes negative and starts to oscillate, so the solution looks 
unphysical. 

If c = 0, it is instead easy to see that eq. (4.8) has a solution for all £, with iV — > as 
£ — > oo. If, on the contrary, c < the real saddle is again lost as iV decreases because of 
the change of sign in the derivative of dz/dN . However, if one chooses Mq ^ | the solution 
eq. (4.3) acquires an extra iV-dependent prefactor, and it can be shown that one recovers 
a stable asymptotic behaviour if M is chosen so that (p(M ) = 0. Hence, also in these 
cases one recovers the generic behaviours found from the Bessel function representation of 
G(£,M)eq. (3.24). 

At first sight the unphysical oscillatory behaviour of G(N, t) when c > appears 
disastrous [21,27-32]. However, here we are not directly interested in the solution eq. (4.3), 
but rather in the dual anomalous dimension which can be extracted from it. This is given 
by 

/ / x d , \ 1 f2/3 N\ 1/3 Ai'[z(aJt), N)} lt . 

7A(a<(t ),i V ) = _i nGWt ) = _ + ^j Ai^CT' (4 ' n) 

which we now study. The large-iV behaviour of 7^4 is easily determined, because when N 
is large, z is large and we can use the asymptotic expansion eq. (4.7) in eq. (4.11). This 
gives 

74a ,, W) .l_ 1 /If^ +0 (_i_). (4 . 12 ) 



2 y k 



Note that this leading behaviour coincides with the 'naive dual' anomalous dimension: 
namely, the anomalous dimension which is found from (p eq. (4.1) using the leading-order 
duality relation (2.9), but with a s = a s (t). We will come back to this point in the next 
section. 

In order to determine the small- N behaviour of we must distinguish three cases. 
For positive values of c, as iV decreases 7 J 4(a s (£), N) increases indefinitely until it blows 
up for the value iV = AT (t) that corresponds to the first zero zq —2.338. . . of Ai(z). 
Because Ai' is regular at z = zq, at iV = Nq the anomalous dimension has a simple pole. 
The location of this rightmost singularity of the anomalous dimension can be determined 
by solving the implicit equation z[a s (t), N ] = zo, with z[a s (t), N Q ] given by eq. (4.5). For 
small values of a s (t) this solution can be expanded as 

iVo(a s , z) = ca s [l + a + 2 -a 2 + V + Q(a% (4.13) 
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Figure 2: Location of the rightmost singularity of the Airy anomalous dimension (4.11), 
No(t) (solid) and of the naive dual anomalous dimension (4.12), a s (t)c (dotted). In 
both cases, the curves correspond (top to bottom) to c = 2, 1.5, 1, 0.5. 

where 

a(a s ,^)^/W /3 (|;) 1/3 - (4-14) 

Recalling that z — > — oo as A" — > and zq < 0, eq. (4.13) implies that Nq < a s c. A plot 
of the position of the rightmost singularity No(t) for several values of c is displayed in 
fig. 2, along with the position of the rightmost singularity of the naive dual anomalous 
dimension eq. (4.12), namely the square-root branch cut at N = a s (t)c. 3 Below Nq, as 
N — > 0, z goes through an infinite number of zeros of the Airy function, and the anomalous 
dimension oscillates wildly. A plot of the anomalous dimension along the positive real axis 
(with c = 1) is displayed in fig. 3: at small A" the anomalous dimension goes through an 
infinity of poles as A" — > 0, while at large N the onset of the square-root drop eq. (4.12) is 
seen. 

If c < 0, it is easy to see that Ai[z(a s (t), N)] does not have any zeros in the complex A" 
plane. Indeed, the zeros of Ai(z) are all located on the negative real z-axis. However, the 
cut induced by the function z(a s (t),N) (4.5) must be taken along the negative A^-axis, so 
that the A^-Mellin inversion integral is well-defined in the physical region N > 0. But the 
function z(a s (t), N) then maps the cut A^ plane on the quadrant — ^ < z < (if c < 0) 
or — ^ < z < ^ (if c = 0). In either case, all the zeros of the Airy function lie outside 
the physical sheet of the A^ plane. Furthermore, in both cases the anomalous dimension 
eq. (4.12) turns out to have a cut along the negative A^ axis starting at A^ = 0, which 
determines the small x asymptotics. 

Specifically, if c = 0, then z oc N 1 / 3 for all N. Because = -3V3E(2/|) the 



3 In this and subsequent figures we take the standard two-loop form of a s (t), with a s (M z ) = 
0.119. 
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Figure 3: The Airy anomalous dimension 7.4(0^, AT) for c = 1 and a s = 0.2 , showing 
the oscillations below the Airy zero at No(t). 



anomalous dimension has a cubic root branch cut at A^ = and is regular for A^ > 0. If 
instead c < 0, the cut can be seen by using the representation of the Airy function in terms 
of the modified Bessel functions K v to rewrite the anomalous dimension (4.12) as 



2(N - a s c)\ 1/2 K 2/3 (jz^) 
^i/3(§^ 2 ) 



It is then easy to show using the continuation formula 

K v (-x) = e~ iU7V K u (x) - ml v (x) 
that t a has a discontinuity along the negative real A" axis, given by (as |AT| 

1/2 



0) 



disc(7A) = — 2i 



-2c 



K 



exp — 



4 1 



3A)|iV|V k 



-2c 3 



+ 0(N). 



(4.15) 



(4.16) 



(4.17) 



Notice that if c < 0, then when A^ is small, z is large and the expansion eq. (4.12) holds. 
This would seem to suggest that the anomalous dimension is regular and positive at A^ = 0. 
However, this conclusion cannot be drawn because this expansion is merely asymptotic: 
in fact, it is factorially divergent, though Borel summable. Its Borel sum coincdes with 
eq. (4.15), with the discontinuity eq. (4.17), which is exponentially subleading and thus 
not immediately visible in the asymptotic expansion eq. (4.12). 

We can now study the splitting function Pa, related by Mellin transform eq. (2.5) to 
the Airy anomalous dimension 7a(ck s (£), AT). If c > 0, the asymptotic £ — > oo behaviour 
of the splitting function is dominated by the pole at A^ = Ao(t), which is the rightmost 
singularity of the anomalous dimension: 



xPA(a s (t),x) ~ a s (t)x~ 



-N (t) 



(4.18) 
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Figure 4: The Airy splitting function Pa(cis,x) with a s = 0.2, and (top to bottom) 
c = 2, 1, 0, —1, —2, showing the absence of oscillations at very small x (solid). The 
standard two-loop splitting function is also shown (dashed). 



If c < 0, the small x behaviour is instead controlled by the cut along the negative 
N axis. Specifically, if c = branch cut is cubic, and the asymptotic behaviour can be 
estimated noting that the inverse Mellin of N 1 ^ 3 is proportional to a gamma function: 
consequently when c = 

xP A (a s (t),x) ~ o ^/3- (4-19) 

If instead c < the asymptotic behaviour can be estimated by evaluating by saddle point 
the integral along the discontinuity eq. (4.17), with the result 

xPA(oi s (t), x) ~ exp 



4 /-2c 3 
V3\ K 



1/4 



1 

A) 



1/2 



(4.20) 



Therefore, asymptotically if c > the splitting function grows as an inverse power 
of x, if c = it drops as a power of £ = In (1/x), and if c < it drops faster than any 
power of £ but slower than a power of x. Whatever the value of c, the splitting function 
is smooth for all £ > 0. A plot of xPA(ot s (t), x) for c = 2, 1, 0, —1, —2 is displayed in 
fig. 4, along with the usual two-loop splitting function: it is apparent that the splitting 
function is well-defined and positive for arbitrarily small x. The fact that the anomalous 
dimension has singularities for N < N (t) (if c > 0) only means that N (t) is the lowest 
moment of the splitting function for which the integral eq. (2.5) over x converges, but 
does not imply that there are problems in the t-evolution. The effect of the running of the 
coupling is to replace the naive dual anomalous dimension eq. (4.12) with the resummed 
anomalous dimension eq. (4.11). We can then equivalently view the effect of the running 
of the coupling as a modification of the evolution kernel, whereby the duality relation (2.9) 
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Figure 5: The kernel (p(M) eq. (4.1) (dashed) and the Airy effective kernel \a{M) 
eq. (4.21) (solid) with c = 1 and a s = 0.2. 

is preserved, but the 'bare' kernel <f(M) (4.1) is replaced by a resummed Airy effective 
kernel, defined by the duality equation 



a s XA(lA(a s , N)) = N 



(4.21) 



in terms of the Airy anomalous dimension 7a(« s , N) (4.11). In the next section we will 
show explicitly that the perturbative effective kernel eq. (2.17) is an (asymptotic) expansion 
of the Airy effective kernel eq. (4.21). The bare and effective kernel (when c = 1) are 
compared in fig. 5. Note that indeed at small M (i.e., by duality, large N) the resummed 
and naive kernel coincide, up to a subleading (asymptotically constant) correction, in 
agreement with eq. (4.12). 

Because the resummed splitting function is smooth, the solution to perturbative evo- 
lution driven by it is also smooth and free of oscillations. This can be understood by noting 
that the solution of the evolution equation (2.3) with the anomalous dimension eq. (4.11) 
and the boundary condition G(N,to) = Fq(N) is 



G(N,t) = T(t,t ;N)F (N), 

W, , M , _ f lA{*s{t),N)dt G(N,t) 

r(Mo;iV) = e ° = GWhV 

with the evolution factor F(t, to; N) explicitly given by 

U ' °' J M[z(a s (t ),N)Y 



(4.22) 



(4.23) 



Now, for large enough ^, the oscillations eq. (4.10) are the same at all values of £, and 
thus cancel in the ratio G(N, t)/G(N, to). This is a consequence of the fact that the saddle 
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Figure 6: Contour plot in the complex iV plane for the Mellin inversion integral of 
G(N,t)/G(N,to). The cross indicates the approximate location of the real saddle, and 
the dotted line is the approximate steepest descent path. 



points which drive the large £ oscillations are dominated by the t-independent —c/N 
contribution to z (4.5), which at small N overwhelms the l/a s (t) term. 

The effect of this i-factorisation is apparent comparing the contour plot of fig. 1 with 
fig. 6, where we display Re [£N + In (G(N, t)/G(N, to))] (all the other parameters being 
as in fig. 1). In this case there is a saddle point on the real axis, while the complex 
saddle points have disappeared. The steepest descent contour goes over the real saddle, 
where the integrand is real, and there are no oscillations. In other words, whenever oscil- 
lations are present they are factorized in the initial condition, while the evolution factor 
r(t,to;N) (4.22) has an inverse Mellin transform which is smooth and asymptotically 
monotonic in £ as £ — > oo. 

In conclusion, in order to understand the impact of the running coupling resummation, 
we give some estimates of the asymptotic behaviour of the inverse iV-Mellin transform of 
the evolution factor F(t,to; N), i.e. of the solution G(N,t) eq. (4.22) when the boundary 
condition Fo(N,t) is constant. If c > 0, the asymptotic behaviour of the evolution fac- 
tor (4.23) is dominated by its rightmost singularity, which is the simple pole at iV = iVo(£o)) 
where the Airy function in the denominator vanishes. The asymptotic behaviour is then 
the same as that of the splitting function eq. (4.19): 

Gfot) ~ aa (t )x- No <- to \ (4.24) 

Note that because No(t) < ot s (t)c the behaviour of the resummed evolution factor is 
significantly softer than that which is found using the naive dual anomalous dimension 
eq. (4.12), namely G ~ x~ aaC (see fig. 2). The same result is of course found by dominating 
the integral with the saddle point displayed in fig. 6. 

If c = the rightmost singularity of the evolution factor is a branch cut at N = 0. In 
this case, the asymptotic behaviour is dominated by a saddle point, located on the positive 
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real axis at 



N„ = 



k(t-t ) 



3/2 



(4.25) 



with k = r(2/3)/r(l/3)3- 2 / 3 (2/3 /K) 1/3 . This leads to 



y/Z 5 / 2 a s (t) 3 / 2 



exp 



k(t - tp) 3 / 2 

£1/2 



(4.26) 



Because the leading asymptotics is constant, in this case we have also given the first 
subleading correction, which corresponds to a power drop in £. This drop is somewhat 
milder than the £~ 3 / 2 drop which is found [23] using the naive dual anomalous dimension 
eq. (4.12). 

Finally, if c < the behaviour can be estimated using a procedure analogous to that 
which led to the corresponding behaviour of the splitting function, eq. (4.20), and in fact 
it turns out to be the same: 



G(£,t) ~ exp 



4 f-2c 3 \ 1/4 n 
V3 V « / \A) 



1/2- 



(4.27) 



Notice that this asymptotic behaviour is subdominant in comparison to that induced by 
a perturbative N = pole, such as is present in the standard one loop anomalous di- 
mension. In fact, if the boundary condition to perturbative evolution has a N = pole 
(corresponding to a constant behaviour in x), the asymptotic behaviour of the evolution 
factor eq. (4.24) for Nq < will be subdominant to that of the boundary condition. It 
is, however, a slower drop in comparison to the power drop G ~ X -Q- S {t)c w hj c ] 1 would be 
found using the naive dual anomalous dimension eq. (4.12) 



5. Asymptotic Behaviour and Resummation: the Generic Case 

The solution of the running coupling ^-evolution equation with a quadratic kernel 
studied in the previous section gives the small x asymptotic behaviour for generic evolution 
kernels x which have a minimum in the physical region < M < 1. Consequently, it is 
possible to use this solution to resum the perturbative running-coupling singularities. 

The key observation which allows one to exploit the results of the Airy resummation 
is that the asymptotic small-a; behaviour of the solution is dominated by the form of the 
kernel in the neighbourhood of its minimum, so we can view the quadratic kernel eq. (4.1) 
as the low-order truncation of the Taylor expansion of the full kernel about its minimum. 
Indeed, it was proven in ref. [23] that if the running is included to finite perturbative order 
(or not at all), then the coefficient of the leading behaviour in the splitting function at 
each order is entirely determined by the value of the kernel and its second derivative at 
the minimum, while higher-order derivatives only affects terms which are asymptotically 
suppressed by inverse powers of £. 
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It is easy to check that this remains true after resummation of the running coupling 
effects (at least if c > 0). To see this, assume that the kernel eq. (4.1) is corrected by a 
cubic term: 



f(M) = c + - K (M - I) 2 + -e(M - I) 3 . (5.1) 



The solution eq. (3.21) can then be written as 



G{N,t) = e~^G (2) (A,t), (5.2) 

where G^(N,t) is the solution which corresponds to a quadratic kernel tp{M), given 
by eq. (4.3). The fact that the cubic term only leads to a subleading correction to the 
asymptotic behaviour is then a direct consequence of the fact that 



d {2f3 N\ 1/3 d 



dt \ k J dz 



(5.3) 



Indeed, assume that the cubic term is treated perturbatively. The change of asymptotic 
behaviour induced by this term is then due to the fact that the position zq of the rightmost 
zero of the Airy solution is shifted by the correction. The shifted zero is located at zq + 5, 
where, expanding eq. (5.3), 5 to first order is given by 

6 K ki{z Q ) V ; 

Hence, 6 is proportional to N and thus suppressed by inverse powers of £ in the 
x — > limit, and so the leading behaviour of the splitting function and solution as x — > 
are still asymptotically given by x~ N °^\ as in the quadratic case eqs. (4. 18), (4. 24). 
A similar line of argument leads to the general conclusion that the leading asymptotic 
behaviour is not affected by corrections to the quadratic approximation, although of course 
subleading corrections to it are affected, as long as c > 0. If c < these arguments fail, 
because the asymptotic behaviour is dominated by the neighbourhood of the point A ~ 
(see eq. (4.20)) i.e. (by duality) the neighbourhood of Mq such %(M) = 0, and not by 
the minimum of x- However in this case running coupling effects are subdominant in 
comparison to the low-order perturbative terms and the boundary condition, (compare 
eq. (4.27)) and thus of no concern. 

We can now exploit the universality of the quadratic result to resum running cou- 
pling singularities, by matching the Airy anomalous dimension to the singular anomalous 
dimension obtained at NLLx. The first step in the procedure is to derive an expansion 
of the Airy anomalous dimension (4.12) of the form eq. (3.20), namely, an expansion in 
powers of a s (t) at fixed a s (t)/N. Because 

[z{a„ N)]- 3 >* = 1 /3 a s 1 , (5.5) 
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the large z asymptotic formula in eq. (4.7) can be used to this purpose. The leading order 
term, already given in eq. (4.12), is now recognized as the leading order of the duality 
expansion (3.20) of the Airy anomalous dimension: 

1 [27N 7 1 p a s 2 
7 ^ s ' N) = 2 - V « U " °\ " 4 1^ + 0( °- ) - (5 - 6) 

So the leading order term has the form that corresponds to ip in eq. (4.1) by the duality 
relation eq. (2.9). The term of order a s (3o in eq. (5.6) is the lowest order effect of running: 
indeed, it is equal to A7 SS eq. (2. 16), (3. 18) for the quadratic ip of eq. (4.1). In fact, 
the asymptotic expansion eq. (4.7) coincides with the asymptotic expansion in powers of 
NPo/ip' 3 eq. (3.19), obtained by evaluating by saddle-point the M-Mellin integral, so the 
expansion of the Airy function eq. (3.20) is recognized to be the same as the running- 
coupling duality expansion of the anomalous dimension eq. (3.20), for the particular case 
of a quadratic x kernel. Hence, the perturbative effective kernel eq. (2.17) is an expansion 
of the Airy effective kernel eq. (4.21). The expansion is factorially divergent but Borel 
summable. 

Furthermore, the universality of the quadratic asymptotic behaviour implies that the 
coefficient of the leading running-coupling singularity to each order is entirely determined 
by knowledge of the value of x an d its second derivative at the minimum. This can be seen 
explicitly at NLLx by considering the running-coupling correction to the dual kernel Axi 
eq. (2.18): whereas of course the function Axi(M) for generic M is determined by the 
full cp, the residue of the pole at M = ^ is entirely determined by knowledge of c, viewing 
eq. (4.1) as the quadratic truncation of the Taylor series of x about M = |: 

Axi = YJT^ + 0{M -^- (5 - 7) 

Therefore, the Airy effective kernel eq. (4.21) is in fact a (Borel) resummation of the leading 
running coupling corections to the kernel eq. (2.17) for generic kernels with a minimum, 
and not only in the particular case of a quadratic kernel. 

It follows that we may resum the running coupling singularities 7 rc eq. (3.20) by adding 
to the anomalous dimension 7 the Airy anomalous dimension which corresponds to the 
quadratic approximation to the kernel x which is dual to 7, and subtracting off the NLLx 
expansion eq. (4.12) in order to avoid double counting. To NLLx it is in fact sufficient to 
determine the Airy resummation from the ^-kernel which is dual to the leading anomalous 
dimension 7 S : explicitly, to this order 



7res (a s , N) = 7s (^) + a s7ss (^) + a s A lss (^ 



+ lA (a Sl N)-^ + ^[N 



a s c + -r 



1 f3 a s N 



(5.8) 



4 N-a s c ' 



where c is the value at its minimum of the x-function dual to j s . Using eq. (5.7) and 
duality, the last term is seen to cancel the singularity in A~f ss (jf), as it must since the 
singularities have been resummed in the Airy anomalous dimension 7^4 eq. (4.11). 
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Figure 7: The leading- (dotted) and next-to-leading order (solid) double-leading x 
kernels compared to the resummed effective kernel (dashed). All curves are computed 
with a s = 0.2 in the R-resummation scheme of ref. [20]. 

The resummed anomalous dimension 7 res eq. (5.8) coincides with the unresummed one 
up to terms which are N 2 LLx, so the resummation does not change the known low-order 
behaviour of the anomalous dimension. In fact, 7^4 is rather small in comparison to the 
standard two-loop anomalous dimension (compare fig. 4), thus it only gives a rather small 
contribution to 7 res (a s , N) eq. (5.8): the main effect of the resummation procedure is the 
subtraction of the running coupling singularities. At asymptotically small x, the leading 
behaviour of the unresummed 7 S and 7 SS , x~ asC , is replaced by the resummed behaviour 
(x~ Ni \ if c > 0). However the subleading unresummed behaviour is not affected by the 
resummation procedure. The subleading behaviour, due to terms in the BFKL kernel 
beyond the quadratic approximation, is suppressed by powers of £ in comparison to the 
leading one, i.e. it is of the form -^x~°' aC . If these terms were resummed, the subleading 
_±_ x -oi 3 c behaviour would be also replaced by a resummed one, x~ N ° +s ^\ with 5 ~ £ -T \ as 
in eq. (5.4). This further subleading running-coupling resummation is lacking in eq. (5.8), 
and it will become necessary when x is so small that the shape of the anomalous dimension 
around N ~ a s c is important. 

By using the resummed anomalous dimension eq. (5.8) in the duality relation eq. (2.9), 
we may finally determine a resummed effective x e ff eq. (2.17). The effect of this resum- 
mation is shown in fig. 7, where we compare the effective x determined at leading and 
next-to-leading order of the double-leading expansion [19,20], and the resummed next- 
to-leading result. The NLO kernel is affected by the singularity, which is then removed 
by the resummation. Note that the resummed result shown in this figure is defined as the 
dual (2.9) of the resummed 7, obtained by combining the Airy resummation eq. (5.8) with 
the double-leading expansion of the anomalous which we discussed in refs. [19,20] (see 
fig. 8 below). For ease of comparison with our previous work, the double-leading kernel 
shown here is the same as that displayed in fig. 2 of ref. [19], which violates the momentum 
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conservation constraint x(0) = 1 by sub-subleading terms. The resummation procedure is 
then applied to this kernel and thus leads to a result which also does not exactly conserve 
momentum. 

In MS-like schemes the running coupling singularities are actually shifted into the 
quark sector, so the resummation should instead be performed in the quark sector. This 
can be done by noting that the running coupling singularities can be removed from the 
anomalous dimension 7+ by redefining the normalization of the gluon distribution by a 
factor u[-f s (a 8 /N)]. Such a redefinition amounts to a conventional change of factorization 
scheme provided only that ^[7(0)] = 1 [33,20]. Upon this redefinition the large eigenvalue 
of the anomalous dimension matrix 7+ (which so far we identified with the anomalous 
dimension, recall eq. (2.1)) and the quark sector anomalous dimension ^ qg change according 
to [33,20] 



+ +' +1/3 Xo(7s) d 
7 + -7 + =7 + + /3o 77 ^- j ^[ln.(M)] 



m= 7s (5.9) 



The choice it(7 s ) = is then sufficient to remove the singular contribution eq. (2.16) 

from 7 + (while the factor in the denominator ensures that this is a bona fide scheme 
change). The scheme change, however, introduces a square-root singularity in the quark 
anomalous dimension r y q9 : indeed, we get 

7 " 7 + H^A?v7~0 V (5 - 10) 

^qgl _ Xoils) ^ qg (5.11) 
\f-Xoils) 



The transformation which takes one to the MS scheme [7] is instead given by a factor 
u(M) = R(M) = r(M)/ s/x^i where r(M) is a regular function of M. Now, we can 
perform the desired running coupling resummation in the quark sector by starting with 
the resummed anomalous dimension eq. (5.8), and then performing a scheme change of the 
form eq. (5.9), while demanding that 7 + be still given by eq. (5.10). The resummation is 
thus shifted entirely into the quark sector, and r ) q9 ' is no longer given by eq. (5.11), but 
rather by 

7 «'(a fl , N) = Ai[z(a s (t), N)]e^ a " N ^~ 3/a ^(a a , N), (5.12) 

which shows explicitly that the square-root singularity in r ) qgl eq. (5.11) is eliminated by 
the resummation (the extra exponential factor removes the term corresponding to 7 S in 
the expansion of the Airy resummed contribution). 

This running-coupling resummation can finally be combined with the small-x resum- 
mation of the anomalous dimension previously discussed by us [19,20]. In such case, the 
expansion of the anomalous dimension in powers of a s at fixed a s /N is obtained by duality 
from a reorganized expansion of the x kernel, where in particular the leading-order term 
is 

a s xo = a s xo + AA, (5.13) 
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Figure 8: The double-leading anomalous dimensions with c = 1 and a s = 0.2 (i.e. 
A = 0.2 [20]) computed to NLO in the R-resummation scheme of ref. [20] with MS 
factorization (dotted). The same after the scheme change eq. (5.10) (solid), displaying 
the running-coupling singularity, and then after the singularity is removed by resum- 
mation (5.8) (dashed). 

with AA = Q!gXi(l/2) + • • • such that A = a s xo(l/2) + AA coincides with the value of 
the all-order x & t its minimum. Consequently, A must be treated as a free parameter. 
To resum the corresponding NLLx running coupling singularity it is sufficient to identify 
the quadratic ip eq. (4.1) with the expansion of xo about its minimum, so that c = X/a s , 
while k is given by eq. (4.2), and then resum according to eq. (5.8). This resummed 
NLLx anomalous dimension can then be combined with the standard two-loop anomalous 
dimension [19,20] to give a fully resummed double-leading expansion. 

The phenomenological impact of these manipulations can be understood on the basis of 
the asymptotic estimates discussed in the previous section. If c > the running coupling 
resummation shifts the rightmost singularity of the small-x anomalous dimension from 
N = a s (t)c to the smaller value N = N (t): N (t) < ot s (i)c, as discussed in the previous 
section (see eq. (4.13)). The asymptotic behaviour is thus somewhat softened. If c = the 
running coupling resummation has a modest effect, in that it simply turns the square-root 
branch point of the unresummed naive dual anomalous dimension eq. (4.12) into a cubic 
branch point; in fact, when c = there are no running coupling singularities, since in this 
case A%i (2.18) and its higher-order generalizations are reg ular asM^| [20]. Finally, 
if c < 0, both the unresummed and resummed behaviours are subdominant in comparison 
to the perturbative behaviour, and thus the phenomenological impact of the resummation 
is negligible. 

This situation is displayed in fig. 8, with the phenomenologically plausible [20] value 
c = 1. In this figure, we show the NLO double-leading anomalous dimension 7+ in the MS 
scheme (used in ref. [19]), and we compare it to the asymptotically unstable anomalous 
dimension obtained from it performing in reverse the scheme change eq. (5.10), so the 



25 



running coupling singularity is moved back into 7+. Note that the momentum conservation 
constraint 7 + (l) = is exactly enforced on all the anomalous dimensions shown in fig. 8, 
by means of a sub-subleading subtraction, as discussed in ref. [19]. The scheme change 
has a very small effect for almost all values of N, but it generates a strong instability in 
a small region of N, 0.2 < N ^> 0.23. (recall that the anomalous dimension has a branch 
cut at N = A) . Finally, we show the result of applying the resummation eq. (5.8) to this 
unstable anomalous dimension: the instability is completely removed, and the anomalous 
dimension remains smooth for all N > X. The softening of the small iV behaviour in the 
vicinity of the branch cut appears to be marginal, as expected because of the smallness of 
the t a contribution to the resummed anomalous dimension eq. (5.8) (see fig. 4). Hence, the 
main effect of the resummation is to remove the instability, thereby widening the region 
of validity of the resummed anomalous dimension down to N = A. Below this value, a 
further subleading running-coupling resummation would be required in order to recover 
the correct asymptotic x~ N ° behaviour. The phenomenological effects of the instability 
and its resummation turn out to be small in the HERA region because of the limited range 
of N in which it appears, but they may become relevant at yet higher energies [20]. 

The running coupling resummation can be generalized to higher orders by using the 
expansion eq. (3.22) in the running-coupling evolution equation, to generate recursive 
equations of the form of eq. (3.23), which determine each Gi in terms of the lower-order 
terms. Note that eq. (3.23) is obtained by neglecting higher order corrections which corre- 
spond to commutators of a s and xo- This terms are nonsingular in the asymptotic limit, 
but will have to be included explicitly into the kernel when going to yet higher orders. It 
is clear from eq. (3.23) that these subleading corrections remain perturbative, in the sense 
that no further singular contributions to x are generated beyond those which appear at 
the leading level. However, the coefficients of these singularities will receive subleading 
corrections which can be determined by explicit solution of eq. (3.23) and its higher order 
generalizations. 



6. Conclusions 

In this paper, we have shown how to introduce the running of the coupling to all 
orders in a BFKL-like x-evolution equation, in a way which is consistent with standard 
renormalization-group evolution. We have thus been able to establish three main results. 

First, we have proven that the solution to the running coupling BFKL equation also 
satisfies a renormalization-group equation, with factorized boundary conditions and a scale 
dependence entirely determined by an anomalous dimension which is related by duality to 
the BFKL kernel. This generalizes to the case of all-order running coupling the duality 
previously derived by us in perturbation theory. Second, by using the well-known Airy 
solution of this running-coupling BFKL equation with a quadratic kernel, we have shown 
that the resummed splitting function is free of the unphysical oscillations which characterize 
the small x behaviour of the solutions to the running-coupling BFKL equation, since these 
get factorized into the initial condition. Finally, we have shown how the Airy solution can 
be used to construct a general resummation of running-coupling singular terms, thanks 
to the fact that the quadratic kernel provides an approximation to the full kernel which 
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correctly determines the asymptotic small x behaviour. The net effect of the resummation 
is small in phenomenologically interesting regions, but it removes a perturbative instability 
which spoiled the asymptotic small x behaviour of the unresummed results. 

In summary, even though it may turn out that in the Regge limit the perturbative 
behaviour is overwhelmed by power-suppressed contributions, or corrections which go 
beyond perturbation theory, we have shown that this need not necessarily be so: a leading- 
twist perturbative approach may provide a consistent description of scaling violations in 
the asymptotic small x limit. 
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